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Abstract. It is shown that explanations of atmospheric neutrino anomaly other than 
— Vt oscillations (e.g. decay, decoherence and v^ — Vt — vkk mixing) can be tested 
at future facilities. Stringent tests of CPT invariance in neutrino oscillations can also 
^ ' be performed. 
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00 I INTRODUCTION 

O 

^ . In this talk I would like to (a) review some non-oscillatory explanations for the 

^ \ atmospheric neutrino data, and how they can be distinguished from the conven- 

or' tional oscillatory explanations in future neutrino experiments; and (b) describe 

p briefly the strong limits that can be placed on CPT violation both from existing 

D . data as well as future experiments especially at neutrino factories. 
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II NEUTRINO DECAY, DECOHERENCE AND 



H : EXTRA DIMENSIONS 

Decay 

If neutrinos do have masses and mixings; then in general, in addition to oscillat- 
ing, the heavier neutrinos will decay to lighter ones via flavor changing processes. 
The only questions are (a) whether the lifetimes are short enough to be interesting 
and (b) what are the dominant decay modes. To be specific, let us assume that 
neutrino masses are at most of order eV[l]. 

For eV neutrinos, the only radiative mode possible is z/j Vj + 7. From the 
existing bounds on neutrino magnetic moments, indirect (but model independent) 
bounds on the decay rates for this mode can be derived: 10~^^s~^, 10~^''s~^ and 
10~^^s~^ for Vt-.v^ and respectively [1]. 
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The decay rate for the invisible three body decay — > Vjuu can be written (for 
rrii ^ rrij) as 



2ri2 5 

1927r3 ^ ' 

The current experimental bound on is of 0(100) [2] (the one loop result in SM is 

< 3.10"^^) and thus V for rrii ^ 0{eV) has to be less than 10"^° s"^ 

The possible existence oi a, 1^ — 0, J — and L — massless particle, x (such as 

a Nambu-Goldstone boson of broken family symmetry) leads to a flavor changing 

decay mode: 

^^Pl+X (2) 

By SU (2) L symmetry decays la —>■ Ip+X occur with the same strength. Current 
bounds on BR{ijl — > ex) and on BR{t nx) of 2.10"^ and 7.10~^[3] respectively 
are sufficient to constrain i/^ and u^- lifetimes to be longer than 10^^ s and 10^° s. 

The only possibility for fast invisible decays of neutrinos seems to lie with ma- 
joron models. There are two classes of models; the 1=1 Gelmini-Roncadelli[4] 
majoron and the 1=0 Chikasige-Mohapatra-Peccei[5] majoron. In general, one can 
choose the majoron to be a mixture of the two; furthermore the coupling can be to 
flavor as well as sterile neutrinos. The effective interaction is of the form: 



Qal^pl^a J (3) 

giving rise to decay: 

J^a^i^p +J (4) 

where J is a massless J = L = 2 particle; Ua and up are mass eigenstates which 
may be mixtures of flavor and sterile neutrinos. Models of this kind which can give 
rise to fast neutrino decays and satisfy the bounds below have been discussed by 
Valle, Joshipura and others [6]. These models are unconstrained by fi and r decays 
which do not arise due to the AL = 2 nature of the coupling. The 1=1 coupling is 
constrained by the bound on the invisible Z width; and requires that the Majoron 
be a mixture of 1=1 and 1=0. The couplings of and i/g {g^ and ge) are constrained 
by the limits on multi-body tt, K decays tt — > nvi'i' and K — > //i/i/i/ and on // — e 
university violation in tt and K decays [7]. 

Granting that models with fast, invisible decays of neutrinos can be constructed, 
can such decay modes be responsible for any observed neutrino anomaly? 

We assume a component of Va, i-e., u^, to be the only unstable state, with a 
rest-frame lifetime tq, and we assume two flavor mixing, for simplicity: 

Ua — cos6i'2 -\- sinOui (5) 

with 7712 > ^1- From Eq. (2) with an unstable 1^2, the Va survival probability is 
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Pcca — sin^O + cos^9exp{—aL/E) 

+ 2sin'^ecos'^eexp{-aL/2E)cos{6m'^L/2E), 



(6) 



where 6m^ = rn^ — m\ and a = 1112 /tq. Since we are attempting to explain neutrino 
data without oscillations there are two appropriate limits of interest. One is when 
the Sw^ is so large that the cosine term averages to 0. Then the survival probabihty 
becomes 



P^^ = sin*e +cos^9exp{-aL/E) (7) 

Let this be called decay scenario A. The other possibility is when 5m^ is so small 
that the cosine term is 1, leading to a survival probability of 

P^^ = {siri^e + cos'^9exp{-aL/2E)f (8) 

We note in passing that scenario A does not provide an acceptable fit to atmo- 
spheric neutrino data [8,9]. Turning to decay scenario B, consider the following 
possibility [10]. The three states f^it, i^r, I's (where Ug is a sterile neutrino) are related 
to the mass eigenstates 1^2, J^s, by the approximate mixing matrix 



/ cos^ sin^ 

-sin^ cos^ I i/3 I (9) 
\ 1 




and the decay is 1^2 — >^ + J ■ The electron neutrino, which we identify with i^i, 
cannot mix very much with the other three because of the more stringent bounds 
on its couplings [7], and thus our preferred solution for solar neutrinos would be 
small angle matter oscillations. 

In this case the 5m\.^ in Eq. (6) is not related to the 5m\j^ in the decay, and 
can be very small, say < 10~^eV^ (to ensure that oscillations play no role in the 
atmospheric neutrinos). Then the oscillating term is 1 and Pi^ij, — > y^) is given by 
Eq. (8). 

The decay model of Equation (8) above gives a very good fit to the Super-K data 
[11] with a minimum = 33.7 (32 d.o.f.) for the choice of parameters 

T^/m^ = 63 km/GeV, cos^ 9 = 0.30 (10) 

and normalization /3 — 1.17. 

The fits (as shown in Fig. 1 in Ref. 10) show the ratios between the Super-K 
data and the Monte Carlo predictions calculated in the absence of oscillations or 
other form of 'new physics' beyond the standard model. The best fits of the two 
models (viz. i/^ — i/^ oscillations and decay) are of comparable quality. The reason 
for the similarity of the results obtained in the two models can be understood 
by looking at Fig. 1, where I show the survival probability Pii^^ v^) of muon 
neutrinos as a function oi L/ E^, for the two models using the best fit parameters. 
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In the case of the neutrino decay model (thick curve) the probabiUty ^(t'^ i'^) 
monotonically decreases from unity to an asymptotic value sin^ 9 ~ 0.49. In the 

case of oscillations the probability has a sinusoidal behaviour in LjEy. The two 
functional forms seem very different; however, taking into account the resolution in 
LjEy^ the two forms are hardly distinguishable. In fact, in the large LjEy region, 
the oscillations are averaged out and the survival probability there can be well 
approximated with 0.5 (for maximal mixing). In the region of small L/E^ both 
probabilities approach unity. In the region LjEy around 400 km/GeV, where the 
probability for the neutrino oscillation model has the first minimum, the two curves 
are most easily distinguishable, at least in principle. 

For the atmospheric neutrinos in Super-K, two kinds of tests have been proposed 
to distinguish between Vy-Vr oscillations and v^J-Vs oscillations. One is based on the 
fact that matter effects are present for v^-Vg oscillations [12] but are nearly absent 
for fy-Vr oscillations [13] leading to differences in the zenith angle distributions 
due to matter effects on upgoing neutrinos [14]. In our case since the mixing is 

— Vj. no matter effect is expected; and hence the recent Siiper-K results [15] are 
in accord with expectations of this decay model. The other test is based on the 
neutral current rate (as measured via production or multi-ring events) which is 
unaffected in i/^ — v-j- oscillations but reduced in — Vs oscillations [16]. In our 
case of the decay model, the neutral current rate is affected and the expectation is 
closer to — i/s mixing. 

Long-Baseline Experiments 

The survival probability of function of L/E is given in Eq. (1). The 

conversion probability into is given by 

P{^^. ^ ^r) = sin2 6 cos^ e{l - e-"^/2^)2 . (11) 

This result differs from 1— P(i/j^ — > u^) and hence is different from v^-Vt oscillations. 
Furthermore, i-*(^'/i — > I'n) + P{i'y — > Vt) is not 1 but is given by 

Pi'^i. ^ J^m) + ^(^M ^^r) = l- cos^ e{l - e-"^/^) (12) 

and determines the amount by which the predicted neutral-current rates are affected 
compared to the no oscillations (or the v^^-Vr oscillations) case. Fig. 2 shows the 
results for ^(t'^ t'^), -P(t'^ Vt) and Piy^j, v^j) + -P(t'^ i't) for the 
decay model and compare them to the v^r'^r oscillations, for both the K2K[16]and 
MINOS [17] (or the corresponding European project [18]) long-baseline experiments, 
with the oscillation and decay parameters as determined in the fits above. 

The K2K experiment, already underway, has a low energy beam E^ ^ 1-2 GeV 
and a baseline L = 250 km. The MINOS experiment will have 3 different beams, 
with average energies Ei, = 3, 6 and 12 GeV and a basehne L — 732 km. The 
approximate L/Ei, ranges are thus 125-250 km/GeV for K2K and 50-250 km/GeV 
for MINOS. The comparisons in Figure 2 show that the energy dependence of z/^ 
survival probability and the neutral current rate can both distinguish between the 
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decay and the oscillation models. ICANOE and especially MONOLITH can also 
test for the oscillation dip[19]. 

Decohcrcncc[20] 

There arc several different possibilities that can give rise to decoherencc of the 
neutrino beam. An obvious one is violation of quantum mechanics, others are 
unknown (flavor speciflc) new interactions with environment etc[21]. Quantum 
gravity effects are also expected to lead to effective decoherence[22,23]. 

The density matrix describing the neutrinos no longer satisfies the usual equation 
of motion: 

p^-i[H,p] (13) 

but rather is modified to 

p=-i[H,p]+D{p) (14) 

Imposing reasonable conditions on L'(p)[24] it was shown by Lisi et al.[20] that 
the Vn survival probability P^^ has the form: 

P^^ = cos^2e + sin'29 e-^'^cos I J . (15) 

where 7 is the decoherence parameter. If is very small {5m?L/2E -C 1), this 
reduces to 

Pn^l = cos^29 + sin^ 29 e'^^ (16) 

If 7 = a/E with a constant, then an excellent fit to the Super-K data can be 
obtained with 9 = 7r/4 and a ~ 7.10"^ GeV/Km. (If gamma is a constant, no fit 
is possible and gamma can be bounded by 10~^^ GeV). The fits to Super-K data 
are shown in ref. 20. and they are as good as the decay or — i/^ oscillations [2 5]. 
The shape of P^^ as a function of L/E is very similar to the decay case as shown 
in Fig. 1. 

Large Extra Dimensions 

Recently the possibility that SM singlets propagate in extra dimensions with 
relatively large radii has received some attention [2 6]. In addition to the graviton, 
right handed neutrino is an obvious candidate to propagate in some extra dimen- 
sions. The smallness of neutrino mass (for a Dirac neutrino) can be linked to this 
property of the right handed singlet neutrino[27]. The implications for neutrino 
masses and oscillations in various scenarios have been discussed extensively [28,29]. 
I focus on one particularly interesting possibility for atmospheric neutrinos raised 
by Barbieri et al [30] . The survival probability P^n is given by 

p,,{L)=\i:l,v.,v:,A{L) r. (17) 

where 
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A{L) = ^ZoUon^'^ exp{iXn^^ L/2ER^) (18) 

where n runs over the tower of Kaluza-Klein states, \^'^\/R^ are the eigenvalues 
of the mass-squared matrix and U'^{^ I/tt^^^) are the matrix elements of the 
diagonalizing unitary matrix. 

An excellent fit to the atmospheric neutrino data can be obtained with the fol- 
lowing choice of parameters: 

^3 = m^R - 3, 1/i? - IQ-^eV, Vl^ ^ 0.4. (19) 

The fit to Super-K data is shown in Ref.27 and obviously it is as good as oscillations. 
This case corresponds to z/^ oscillating into z/^ and a large mimber(about 25) of 
Kaluza-Klein states. Because of the mixing with a large number of closely spaced 
states, the dip in oscillations gets washed out and P^^ looks very much like the 
decay model as shown in the Fig. 1. 



Ill CPT VIOLATION IN 
NEUTRINO OSCILLATIONS [31] 

Consequences of CP, T and CPT violation for neutrino oscillations have been 
written down before [32]. We summarize them briefly for the Ua — > i^p flavor 
oscillation probabilities Pq^ at a distance L from the source. If 

P^p{L) ^ P^-^iL) , P^a, (20) 
then CP is not conserved. If 

Pap{L) ^ Pf,^{L) , P^a, (21) 
then T-invariance is violated. If 

P^p{L) ^ P-p^{L) , P^a, (22) 

or 

Paa{L) ^ Paa{L) , (23) 

then CPT is violated. When neutrinos propagate in matter, matter effects give 
rise to apparent CP and CPT violation even if the mass matrix is CP conserving. 

The CPT violating terms can be Lorentz-invariance violating (LV) or Lorentz 
invariant. The Lorentz-invariance violating, CPT violating case has been discussed 
by CoUaday and Kostelecky [33] and by Coleman and Glashow [34] . 

The effective LV CPT violating interaction for neutrinos is of the form 

^2bflA , (24) 
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where a and (3 are flavor indices. We assume rotational invariance in the "preferred" 
frame, in which the cosmic microwave background radiation is isotropic (following 
Coleman and Glashow [34]). 



(25) 



where is a hermitian matrix, hereafter labeled h. 

In the two-flavor case the neutrino phases may be chosen such that h is real, in 
which case the interaction in Eq. (24) is CPT odd. The survival probabilities for 
flavors a and a produced at t = are given by [34] 



and 
where 



Paa{L) = 1 - sin^ 29 sin2(AL/4) , 

Paa(L) = 1 - sin2 2esin2(AL/4) , 

A sin 20 = I (5m^/E) sin 2^^ + 25be''' sin 2^^ 
A cos 26 = (Sm'^/E) cos 29^ + 25b cos 2^^ ■ 



(26) 

(27) 

(28) 
(29) 



A and are defined by similar equations with 6b — > —6b. Here 9^ and 9b define 
the rotation angles that diagonalize and b, respectively, = m| — and 
6b — b2 — bi, where mj and bi are the respective eigenvalues. Wc use the convention 
that cos2^^m and cos 2^5 are positive and that 6m'^ and 6b can have either sign. 
The phase f] in Eq. (28) is the difference of the phases in the unitary matrices 
that diagonahze and 6b; only one of these two phases can be absorbed by a 
redefinition of the neutrino states. 

Observable CPT- violation in the two-flavor case is a consequence of the interfer- 
ence of the terms (which are CPT -even) and the LV terms in Eq. (24) (which 
are CPT-odd); if 6w? = or 56 = 0, then there is no observable CPT- violating 
effect in neutrino oscillations. If 6w?/E » 26b then © ~ and A ~ 6vr? jE., 
whereas if 6w? jE -C 26b then ~ and A ~ 26b. Hence the effective mixing an- 
gle and oscillation wavelength can vary dramatically with E for appropriate values 
of 6b. 

We note that a CPT-odd resonance for neutrinos (sin^ 20 = 1) occurs whenever 
cos 20 = or 



{6m^/E) cos 29m + 26b cos 29b = 0; (30) 

similar to the resonance due to matter effects [35,36]. The condition for antineutri- 
nos is the same except 6b is replaced by —6b. The resonance occurs for neutrinos 
if (5m^ and 6b have the opposite sign, and for antineutrinos if they have the same 
sign. A resonance can occur even when 9^, and 9^ are both small, and for all values 
of T); if 9m — 9b, a resonance can occur only ii rj ^ 0. If one of ly^ or i/^ is 1/^, then 
matter effects have to be included. 
If 77 = 0, then 
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A= {6m'^/E) + 25b. 



(31) 
(32) 



In this case a resonance is not possible. The oscillation probabilities become 



For fixed E, the 5b terms act as a phase shift in the oscillation argument; for fixed 
L, the 5b terms act as a modification of the oscillation wavelength. 

An approximate direct limit on 5b when a = fi can be obtained by noting that in 
atmospheric neutrino data the flux of downward going z/^ is not depleted whereas 
that of upward going t/^ is[ll]. Hence, the oscillation arguments in Eqs. (33) and 
(34) cannot have fully developed for downward neutrinos. Taking \5bL/2\ < 7r/2 
with L ~ 20 km for downward events leads to the upper bound \5b\ < 3 x 10~^° GeV; 
upward going events could in principle test \5b\ as low as 5 x 10^^^ GeV. Since the 
CPT-odd oscillation argument depends on L and the ordinary oscillation argument 
on L/E, improved direct limits could be obtained by a dedicated study of the energy 
and zenith angle dependence of the atmospheric neutrino data. 

The difference between P„„ and Prt^i 



can be used to test for CPT- violation. In a neutrino factory, the ratio of z/^ z/^ 
to — > z/^ events will differ from the standard model (or any local quantum field 
theory model) value if CPT is violated. Fig. 3 shows the event ratios Ni^u^ — > 
z/^)/iV(z/^ u^) versus 56 for a neutrino factory with 10^^ stored muons and a 10 kt 
detector at several values of stored muon energy, assuming (5m^ = 3.5 x 10^'^ cV^ and 
sin^2^ = 1.0, as indicated by the atmospheric neutrino data[ll]. The error bars in 
Fig. 3 are representative statistical uncertainties. The node near 5b = 8x10"^^ GeV 
is a consequence of the fact that Paa — Paa, independent of E, whenever 5bL — mr, 
where n is any integer; the node in Fig. 3 is for n = 1. A 3a CPT violation effect 
is possible in such an experiment for 5b as low as 3 x 10~^^ GeV for stored muon 
energies of 20 GeV. Although matter effects also induce an apparent CPT- violating 
effect, the dominant oscillation here is z/^ — > i>t, which has no matter corrections 
in the two-neutrino limit; in any event, the matter effect is in general small for 
distances much shorter than the Earth's radius. 

We have also checked the observability of CPT violation at other distances, 
assuming the same neutrino factory parameters used above. For L = 250 km, the 
5bL oscillation argument in Eq. (35) has not fully developed and the ratio of P to u 
events is still relatively close to the standard model value. For L — 2900 km, a 5b 
as low as 10~^^ GeV may be observable at the Scr level. However, longer distances 
may also have matter effects that simulate CPT violation. 




(33) 



(34) 




(35) 
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IV SUMMARY 



At Long Baseline Experiments and Neutrino Factories true signatures of oscil- 
lations (dips) can be established and decay like scenarios can be excluded with 
confidence. Furthermore these facilities can test CPT conservation at levels better 
than 10^3 GeV. 
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FIGURE 1. Survival probality for versus logiQ (L/E) for the decay model, decoherence, extra 
dimensions and oscillation. 
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FIGURE 2. Long-baseline expectations for the K2K and MINOS long-baseline experiments 
from the decay model and the v^-Vr oscillation model. The upper panel gives the neutral current 
predictions compared to no oscillations (or v^-Vr oscillations). 
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FIGURE 3. The ratio of 9^ to event rates in a 10 kt detector for a neutrino 

factory with 10^^ stored muon with energies E'^ — 10, 20, 30, 50 GeV for basehne L = 732 km 
versus the CPT-odd parameter Sb with 9^ ~ db = S and phase 77 = 0. The neutrino mass and 
mixing parameters are Sm^ ~ 3.5 x 10^^ eV^ and sin^ 26* — 1.0. The dotted hne indicates the 
result for 8h = 0, which is given by the ratio of the v and v charge-current cross sections. The 
error bars are representative statistical uncertainties. 
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